Expanding Universe : Thermodynamical Aspects From Different 
Models 



Sridip Pal *0 
Ritabrata Biswas H 



* Department of Physical Sciences, Indian Institute of Science Education and Research, Mohanpur - 74-1252, 
Nadia, West Bengal, India 

f Department of Physics, Indian Institute of Science, Bangalore-560012, India. 

Abstract 

The pivotal point of the paper is to discuss the behavior of temperature, pressure, energy 
density as a function of volume along with determination of caloric EoS from following two model: 
w(z) = w + wi ln(l + z) & w(z) = — 1 + a +2^^i+z)+aI\i+z) 2 ■ T ne tmie sca l e 01 instability 
for this two models is discussed. In the paper we then generalize our result and arrive at general 
expression for energy density irrespective of the model. The thermodynamical stability for both 
of the model and the general case is discussed from this viewpoint. We also arrive at a condition 
on the limiting behavior of thermodynamic parameter to validate the third law of thermodynamics 
and interpret the general mathematical expression of integration constant U$ (what we get while 
integrating energy conservation equation) physically relating it to number of micro states. The 
constraint on the allowed values of the parameters of the models is discussed which ascertains 
stability of universe. The validity of thermodynamical laws within apparent and event horizon is 
discussed. 
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1 Introduction 

To explain the cosmic acceleration predicted from the la type supernova observations (Perlmutter, S. et al. 1999; 
Riess, A. G. et al. 1998) one popular wayout is to modify the stress energy tensor part, i.e., the right hand side 
of the Einstein's field equation. Existence of some unknown matter termed as Dark Energy(DE hereafter) is 
been assumed (Riess, A. G. et al 2004; Perlmutter, S. et al 1998; Garnavich, P. M. et al 1998; Bachall, N. A. 
et. al. 1999; Copeland, E. J. et. al. 2006) which violates the strong energy condition. The simplest candidate 
of DE is a tiny positive cosmological constant (A) which obeys the equation of state (EoS hereafter), w = — 1. 
But due to low energy scale than the normal scale for constant A, the dynamical A was introduced (Caldwell, 
R. R. et. al. 1998). Again at very early stage of universe the energy scale for varying A is not sufficient. So to 
avoid this problem, known as cosmic coincidence (Steinhardt, P. J. et. al. 1999), a new field, called tracker field 
(Zlatev, I. et. al. 1999) was prescribed. In similar way there are many models (Sahni, V. and Starobinsky, A. A. 
2003) in Einstein gravity to best fit the data. Yet its require some modifications. From this point of view some 
alternative models are evolved. Most of the DE models involve one or more scalar fields with various actions 
and with or without a scalar field potential (Maor, I. et. al. 1998). Now, as the observational data permits us 
to have a rather time varying EoS, there are a bunch of models characterized by different scalar fields such as 
a slowly rolling scalar field (Quintessence) (—1 < w < —1/3, w(= p/p), being the EoS parameter) (Caldwell, 
R. R. et. al. 1998), k-essence (Armendariz - Picon et. al. 2000), tachyon (Sen, A. 2002), phantom (uj < — 1) 
(Caldwell, R. R. 2002), ghost condensate (Arkani-Hamed, N. et. al. 2004; Piazza, F and Tsujikawa, S. 2004), 
quintom ( Feng, B. et. al. 2005), Chaplygin gas models (Kamenshchik, A. Y. et. al. 2001) etc. Some recent 
reviews on DE models are described in the ref. (Copeland, E.J. et. al. 2006; Li, M. et. al.2011). 

While explaining evolution of the universe, various DE models have been proposed, all of which must be 
constrained by astronomical observations. In all the models, the EoS parameter ui plays a key role and can 
reveal the nature of DE which accelerates the universe. Different equations of state lead to different dynamical 
changes and may influence the evolution of the universe. The EoS parameter uj and its time derivative with 
respect to Hubble time are currently constrained by the distance measurements of the type la supernova and 
the current observational data constrain the range of equation of state as — 1.38 < to < —0.82 (Melchiorri, A. 
et. al. 2003). Recently, the combination of WMAP and Supernova Legacy Survey data shows a significant 
constraint on the EoS uj = — 0.971q]q9 f° r the DE, in a flat universe (Seljak, U. et. al. 2006). Recently some 
parametrization for the variation of EoS parameters ui(z) have been proposed describing the DE component. 
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1. lj(z) = cj + co\z (Cooray, A. R. and Hutcrcr, D. 1999). Here w = — 1/3 and wi = —0.9 with z < 1. 
This grows increasingly unsuitable for z > 1. So the following model has been proposed. We will call this 
parametrization as 'linear parametrization' . 

2. oj(z) = ljq + wiy^-j. This ansatz was first discussed by Chevallier and Polarski (Chevallier, M. and 
Polarski, D. 2001) and later studied more elaborately by Linder (Linder, E. V. 2003). The best fit values 
for this model while fitting with the SNIa gold dataset are uj = —1.58 and wi = 3.29. We will call this 
parametrization as l CPL{Chevallier — Polarski — Linder) parametrization' . 

3. oj(z) = u>o + u>i (i+z)2 ( Jassal, H. K. et. al. 2005). A fairly rapid evolution of this EoS allowed so that 
uj(z) > — 1/2 at z > 0.5 is consistent with the supernovae observation. We will call this parametrization as 
l JBP{Jassal — Bagla — Padmanabhan) parametrization' . 

4. lj(z) = -1 + A +2A 1 1 a+t 2 )+Xli+z) 2 ( Alam > U - et - aL 2004a, 2004b). This ansatz is exactly 
the cosmological constant lj = — 1 for A\ = A 2 = and DE models with w = —2/3 for Ao = A2 = 
and to = —1/3 for Aq = A\ = 0. It has also been found to give excellent results for DE models in which 
the equation of state varies with time including quintessence, Chaplygin gas, etc. The best fit values of A\ 
and A 2 are A\ = —4.16 and A 2 = 1.67 for the SNIa Gold dataset. We will call this parametrization as 
i ASSS(Alam — Sahni — Saini — Starobinski) parametrization' . 

5. oj(z) = cjo + ln(l + z) (Efstathiou, G. 1999). This evolution form of EoS is valid for z < 4. lji is a 
small number which can be determined by the observations. The minimum value of wi is approximately —0.14 
and wo > — 1- We will call it 'Log parametrization' . 

Thermodynamics of DE universe filled up with the fluids with linear and JBP parametrization had been 
studied before (Mazumder, N. et. al. 2010). Thermodynamics with CPL was studied by Xing, L. et. al. (2011). 
In this paper we are going to study DE universe filled with fluid Log parametrization and ASSS parametrization 
respectively which have not been studied before so extensively. 

As we went through the literature we came to know that apart from the cosmological constant (w = —1), the 
sinchypcrbolic scalar field potential (Sahni, V. and Starobinsky, A. A 2000; Urena-Lopez, L.A. and Matos, T. 
2000; Sahni, V. et al. 2003) and the topological defect models, there is no DE model with constant w consistent 
with observation (Alam, U. et. al. 2004a). Almost every model like Quintessence, Chaplygin gas (Sctarc, M.R. 
and Bouhmadi-Lopez, M.R. 2007, Setare, M.R. 2007b) etc. depicts significant evolution in uj(z) over sufficient 
time scale. 

The way towards a meaningful reconstruction of tj(z) depends on inventing an efficient fitting function for 
either d[ y (z) or H(z), where c?l(z) is luminosity distance and H (z) is Hubble parameter. The parameters of this 
fitting function are determined by matching to Supernova observations. Now, one can manage to reconstruct 
the functional form of u(z) by taking an ansatz of the luminosity distance or H{z) and then by comparing it 
to supernova observations. 

This model is based on the following ansatz: 



where x = (1 + z), po c is the critical density at present time and f2 m dimensionless density of dark energy. 

If we carefully note we can see the ansatz in terms of energy density is actually a truncated taylor expansion 
where it is exploited that any well bahaved function can be well approximated by Taylor expansion within a 
range. So the model is valid for z < few (Alam U et. al. 2004a). 

The model yields excellent results among DE models in which the EoS varies with time including quintessence, 
Chaplygin gas etc. (Sahni, V. et al. 2003; Alam, U. et al. 2003). Alam, U et al. 2004a (In figure 1 of their paper) 
have shown the accuracy of the ansatz when it is applied to several other DE models like tracker quintessence, 
the Chaplygin gas and super-gravity (SUGRA) models. They also plotted the deviation of log (diiJo) (which 
is the measured quantity for SNe) obtained via the ansatz from the actual model values. Clearly the ansatz 
comes out to be of excellent agreement well over a significant red-shift range for f2o m =0.3. The ansatz agrees 




(1) 



Which can be equivalently written as following: 



Pde = Poc (A + A x x + A 2 x 2 + A 3 x 3 ) 



(2) 
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with these models of dark energy with less than 0.5% errors in the redshift range where we do have SNe data 
available (Alam, U. et. al. 2004a) 

The analysis of type la supernova data involving the priors with most frequently used condition as uj = 
constant and uj > 1 leads to confinement of DE to within a narrow class of models. Moreover, when we impose 
such priors on the cosmic equation of state it can result to a complete misrepresentations of reality as shown in 
(Maor, et al. 2002). 

Recently, evolution of EoS parameter with red shifts been investigated (Vazquez, J. A. et. al. 2012) by 
performing a Bayesian analysis of current cosmological observations. Vazquez, J. A. et. al. have shown if they 
calculate the Bayes' factor then most of the data catalogues supports the fact that CPL, JBP are significantly 
disfavored with respect to simple AC DM model (uj = —1). Now, linear parametrization is not to be considered 
in such a case as the parametrization blows with high z. Calculation of Bayes' factor for ASSS and kCDM 
using codes like CAMB etc are bit cumbersome. But from intuition we can see —1 term is already there in 
ASSS parametrization. Extra additional z dependent term may make it different enough from the ACDM 
model. Even theoretically, ASSS is efficient enough to explain different phases of universe. These observational 
supports substantiate our interest of working with ASSS model. 

Linear, CPL, JBP or ASSS - all the parametrizations can be treated as kind of combinations of polynomials 
with (1 + z) as polynomial parameter. These modes are efficient for early epochs whereas they are not that 
effective for low redshifts. While parametrizing, the fact that at low redshifts the magnitude-redshift relation 
is degenrate for the models having same deceleration parameter, should be taken into account. Besides keeping 
consistency of SN constraints on dimensionless DE density with those derived from CAMB measurements is 
another important factor. These constrains along with many others had given birth of Log parametrization. 
This is somehow not following the common pathway of considering DE EoS. These properties did motivate us 
to take Log parametrization for the study of thermodynamics of universe. 

We will study the universe from thermodynamical aspect for both type of fluids separately in Sec. [2] In 
Sec. [3] we will seek the thermal EoS for both type of fluid separately along with study on thermodynamical 
stability. The Internal Energy, thermodynamic pressure, entropy, temperaure as a function of volume have 
been determined in this section. The Sec. [4] contains a general derivation of thermal EoS for any fluid obeying 
p = uj(z)p and a general discussion on stability criterion. We have tried to give a physical interpretation of 
the generality of thermal EoS and the expression for energy density as a function of temperature and volume 
in this section. Study of instable cases, onset of instability have been discussed in Sec. [5] The Sec. [6] proves 
the validity of laws of thermodynamics on the apparent horizon and invalidity of the same on the event horizon 
which agrees with (Mazumder, N. et. al. 2010; Xing, L. et. al. 2011). The paper ends with a brief concluding 
remark in Sec. [7] 

2 Study of Universe Treating it as an Thermodynamical System 

Let us consider an Universe filled with a perfect fluid having volume V and p, p, T and S are respectively 
the energy density, thermodynamical pressure, temperature and entropy of the system. From the first law of 
thermodynamics (Myung, Y. S. 2011) we have 



TdS = d(pV) + P dV = d[(p +p)V]- Vdp 



(3) 



We have the following integrability condition (Myung, Y. S. 2011), i.e., 



d 2 S d 2 S 



dTdV dVdT 



which yields (Gong, Y. et. al. 2007) 



dp dT 



(4) 



P + P T ' 



So combining equations Q and Q and integrating we obtain the following: 



S = 



(p + p)V 



(5) 



T 
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where we have dropped an additive constant devoid of physical significance. 

We assume our universe to be homogeneous and isotropic FRW space-time with following line element 



ds 2 = -dt 2 +a 2 (t) 



r 2 (dO 2 + sin 2 9d6 2 ) 



1 — kr 2 



(6) 



where k, the curvature scalar having values 0, ±1 for flat, closed and open universe respectively. The Friedmann 
equations and the energy conservation equation are 

2 8nGp k 

H = ~3~ ~ 2 ' (?) 

H = -4irG(p + p) + ^ (8) 
or 

and 

P + 3H(p + p) = , (9) 
where the Hubble parameter is given by H = ^ . 

2.1 Thermodynamics of Fluid with u)(z) = Uq + uj\ In (1 + z) 

Using this u)(z), integrating the energy conservation equation (JSJ we would get the expression of energy density 
as a function of redshift as follows : 

p = po(l + z fi l +^+^ l < 1 + z )} , (10) 
where po is integration constant denoting the present time (z = 0) density of our universe. 

1. From the graph of energy density vs redshift, we observe that p initially has a very high value. With 
expansion of the universe, u decreases, hence p falls to po as z becomes from 2. 

2. The graph of Pressure vs Redshift reveals that the pressure remains positive for allowed values of uio and 
U3% in the range z € (0, 2). It initially starts with high value and then decreases down to present value as 
z goes from 2 to 0. 

The integration of the integrability condition ((4]) gives expression for temperature as 

T = T (l + W )(l + z )3{"o + ^/n(l + 2 )} j 

where To is the integration constant (Note To (I + loq) is the present time temperature). 

• The graph between Temperature and Redshift reveals that Temperature is also dropping as we traverse 
from past to present, i.e, as z varies from 2 to 0. 

Using the last two expressions and plugging in the expression for entropy in the last section, we get the following 
expression for S : 

S=£. (12) 
The heat capacity and square of the sound velocity are given by: 

dp + 

c v{z) = V— = — — . 13 

Ol U\ + 6UJ (1 + LI) 

Note when u < — 1 we have cy < i.e universe is in a unstable phase. Hence the stability condition demands: 
uj > -1. 
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Fig 1 

Energy Density vs redshift 




Fig 2 

Pressure vs redshift 




Fig 3 

Temperature vs redshift 




0.0 0.5 



The expression for the sound's speed 



2/ \ dp 



dp 3(1 + u) 



Now demanding u s 2 < 1 we get the following condition: 



{lo + wi/n(l + z)} < 



2 „ 3 - Wi 



3 ' 



along with the constraint w\ < 3. z has a minimum value of 0, hence we can write: 

3wq + oj\ < 3. 



(14) 

(15) 
(16) 



2.2 Thermodynamics of Fluid with W (z) = -1 + 

Using this w(z), integrating the energy conservation equation we would get the following: 

, f / M ( A Q +A 2 (l + z) \] 

p = p Q M exp <^ arctan — V , (17) 

{ ^A A 2 - A\ \ ^A A 2 - A{ J J 

where po is integration constant which we can get putting z = 0, and M = Ao + 2j4i(1 + z) + A 2 (l + z) 2 . 
On the other hand integrating the integrability condition, we get the expression for temperature: 

T = T (l + W )(l + zr^MexJ- f arctan f ^+^+A ) ), (18) 

[ y/A A 2 - A\ \ yjA a A 2 - A\ J J 

where To is the integration constant. 

Using the last two expressions and plugging in the expression for entropy in the last section, we get the 
following expression for S: 

S=£. (19) 
The heat capacity and square of the sound velocity are given by: 

Q ( 1 3 l^i 1 

Upon simplification we get: 

3S(2 + c) , (or , 

cy = n 5— and 21 

3(l+ W ) 2 + (l + z)f 

2 f , dp dw \ A 1+ A 2 (l + z) y 1 

vAz)= d-p =UJ + ^{ A + 2Ml + z) + A 2 (l + zr j (22) 
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3 Derivation of Thermal EoS &; Study of Stability 



In this approach it is very convinient to consider the fluid obeying adiabatic EoS p = u)(z)p with constant 
particle number TV as a thermodynamical system. Without any loss of generality we can assume the internal 
energy (U) and the pressure (p) as functions of entropy (5*) and volume (V). So we can structurize our density 
and pressure and also the concerend differential equation as (Landau, L. D. and Lifschitz, E. M. 1984) 



p = -(w)s and 
=0. 



(23) 



3.1 Fluid with u(z) = lo + u\ In (1 + z) 

In (|2U)) . use of w(z) = uj + u)iln(l + z), yields the following: 

U = U Q V^ lnV) exp (-uoV), (24) 

where Uq is integration constant (It can be function of S and uio should be greater than for U not to diverge). 
Hence the energy density becomes 

p = U V i ^ lnV - 1) e^{-Lu V), (25) 
while the expression for pressure is following: 

p = Uq [uj - Y lnV } V^ lnV ~V exp {-uoV). (26) 
The criterion for stability of the fluid during expansion: 



1. (ft) 0. 



s 

2. The thermal capacity at constant volume should be greater than zero i.e cy > 0. 
The first condition leads to the following constraint: 

(^.I.^V) ^)<^. (27) 



Now to get to the thermal equation of state we start with the expression for temperature: 

dU 
~dS 



Using the expression for internal energy we obtain: 

r =^f = fr»-)^(-^. (29) 

Now we will look into Uq by considering the change of variable from (S,V) to (P,T). The Jacobian of the 
transformation is (Landau, L. D. and Lifschitz, E. M. 1984): 

'-(&).(£).-(*).(£)„• 



Note, we can also write it out as the following: 

dSJ v \ dS 
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Equating the last two expressions and using the expression for U and T we get the following: 



(wo-ylnV) (- Wo + ^lnv) lu Q 



TT d 2 U Q fdU 



dS 2 



dUy 

dS 



= 







=> Uo = cxp (aS), 

where a is integration constant upon integrating the differential equation for Uo- Hence we find 

T 

p=— & 
av 

p = up= — — (cjo - \uiln~V 
av \ 3 

The existence of critical point requires the following condition to be true: 

(*),- 1 
(^A -o 

(^) T <0 



(32) 

(33) 
(34) 

(35) 
(36) 



(37) 



But in this model, the point where we have = ^ nas ^ ne P ro P cr ty that second of above criterion does 

not hold as 

(d 2 P \ _Tu J "(V) _ lo\ T {m 

where primc(/) denotes derivative with respect to volume. Hence, in this model universe does not go through 
any critical point, even if there is a transition from unstable to stable configuration or vice- versa the transition 
is smooth. 

The specific heat at constant volume comes out to be : 



dT J 



(39) 



v 



Choosing a > guarantees the fulfillment of second criterion for stability, i.e., cy > 0. while the velocity of 
sound is given by: 

„.2 



wi 

uiilnV — 3 — 3^0^ 

Here also v 2 < 1 imposes constraint on the allowed values of w and w i • 
From the expression for Uq, the entropy is given by: 



5=1 

a 



3.2 Fluid with u(z) 



1 



1 + 2 



in | -) +uj V - ^{InVf 



A- L +2A 2 (l+z) 



(40) 



(41) 



3 A +2A 1 (l+z)+A 2 (l+z) 2 
The expression for U in this case becomes: 



U = U (a V + 2A 1 V 2/3 + A 2 V 1/3 ^j exp < 



' Al arctan ( A f 1/3+A i) 



^A A 2 - A 2 



(42) 
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where Uq is integration constant. Therefore energy density comes out to be: 



p = U Uo + 2A 1 V- 1 ' 3 + A 2 V- 2 /^ exp 



Ax arctan I A f ' +A \ 



^/A A 2 - A\ 



while the pressure is given by 

V = u](V)U (a + 2A X V- 1/3 + A 2 V- 2/3 ^j exp < 



A\ arctan 



A„V 1 ' 3 +A 1 
y/A A 2 -Af 



y / A A 2 - A\ 



where «(*) = (-1 + and V = a* = (1 + *)"». 

Now we have two criterions for stability of the fluid during expansion 



1. 



(*), 



< 



2. The thermal capacity at constant volume should be nonnegative i.e cy > 
The first condition yields: 



.1.. + 2.i,r - i i \) | f^ + <*>(K) 



where i? = (AxV^ + A 2 V~^ u(V). 

Now we start with the expression for temperature to obtain thermal EoS: 



T 



dU 
~dS 



Using equation (|42|) we obtain: 

U dU dU 



A\ arctan 



T = 



= ^( A oV + ZA.V 2 ' 3 + A.V 1 ' 3 ) exp ■ 



AgV 1 ' 3 +A 1 

y/A A 2 -A\ 



Uq dS dS 



Equating (|30f and (|3Tj) for this system we have, 



V / A A 2 - A\ 



U, 



d 2 U fdU^ 1 



dS 2 \ dS 



= ^ Uq = exp (aS) 



(43) 



(44) 



(45) 



(46) 



(47) 



This result resembles with Even we can see the expression for energy density will be similar to (|35l) . 

Expression for pressure in this case is : 



V = up = < -1 



l + z A 1 + 2A 2 (l + z) \ T 

3 A + 2Ai(l + z) + A 2 {1 + z) 2 J aV 



(48) 



The same conditions (|37j) applies here too which yields the following conditions needed to be satisfied for having 
critical point: 

u"(V) = & u'"(Y) < (49) 
The thermal capacity of the system at constant volume becomes : 



v ,'dS\ 1 
1 = T[ W)v = a 



(50) 
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If we choose a > it automatically guarantees the fulfillment of second criterion for stability in the same way 
as previous section i.e cv > 0. 

While the velocity of sound is given by: 



dV 

Now v 2 < 1 imposes constraint on uig and uj\. 

From the expression for Uq, we obtain the expression for entropy which is: 



S = 



h, | - ) - In ( AoV + 2A 1 V 2/3 + A 2 V 1/3 
a 



Ai arctan 



a v 1 ' 3 +a 1 

y /A A 2 -A 2 1 



VAqA 2 



A\ 



(51) 



(52) 



4 A General Derivation : Independent of Model 

Looking deep into the fact that in last section we have observed Uq comes out to be exp (ccS) independent of 
the form of u)(z), here we are giving a proof without assuming any specific form of u)(z) which makes the result 
stronger. 

We start with integrating energy conservation equation (|23l) 

U = U F(V) (53) 
where F(V) = exp j— J^^p-dV^. Hence P becomes: 

F(V)u){V) 



U - 



V 



Now to get to the thermal equation of state we start with the expression for temperature 
expression for internal energy we obtain: 



T 



U dU dU 



F(V) 



U dS dS 

Equating (f30|) and ([3"Tjl and using the expression for U and T we get the following: 



F(V)F'(V)lu(V) 
V 



u 



d 2 U 
dS 2 



dU 
dS 







(54) 
and using the 

(55) 
(56) 



where a is integration constant upon integrating the differential equation for Uq- Immediately, the expression 
for the energy density becomes p = while the pressure is given by 



_ ujT 
P ~aV 

The specific heat at constant volume comes out to be cy = T (fy) y = 
and the entropy is given by following expression: 



(57) 



5=1 

a 



I - ) lnF(V) 



(58) 



So, T should be > aF(V) to ascertain S > 0. 
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The third law of thermodynamics demands the entropy to go to as temperature approaches zero. This 
demand translates into following criterion when we use (|58|) : 



H— = a =Cv (59) 

So we can conclude F(V) must tend to as universe cools down to absolute zero. Using this information, from 
equation ([59]) we can arrive at the fact that F(V) is actually a measure of change of internal energy to effect a 
small change in temperature around T = point. 

Physically Interpreting The General Form of Uq : 

Irrespective of the model we have obtained, 

U = cxp (aS) (60) 

In its differential form we can write 

— — = adS (61) 

Now note we have S = kslnfl where is the number of micro states corresponding to a macrostate of the 
system. In its differential form it looks very much like (|6Tj) 

So comparing equation (1611) and (|62[) . we interpret Uo as a measure of number of micro states. Now Uq is related 
to energy of universe at present and it is physically very plausible to have energy going in proportional to the 
number of micro states. 



5 Instability : Time scale for Onset of Instabilities 

He was Einstein in very known history who has made "the blunder" while producing a static universe which 
was unstable and not observationally supported! Latter instability regarding cosmological, matter creation has 
been studied (Saslaw, W. C. 1967). General idea of thermodynamics says that negative specific heat indicates 
an instability (of course thermal). 

Logarithmic Model 

Here the expression of cy illustrates two cases of instability: 

(1) to < — 1 along with u)\ + 3w(l + uj) > : Now it can be easily noted that the former criterion implies the 
later. Hence we have instability for uj < —1. Taking ujq = —0.995 and u)\ = 0.25 (Efstathiou, G. 1999) we 
would get that z — —0.019 marks the onset of instability. As z goes in the opposite direction of time, negative 
z indicates towards the future. So z — —0.019 can be treated as the timescale for onset of instability. This may 
suggest that our universe is now in a thermally stable equilibrium but is heading towards an instability where 
small perturbation can lead to catastrophic change. 

(2) — 1 < to along with to% + 3a; (1 + uj) < : Now u > does not satisfy this condition, hence we have 
stability when uj > 0. So the range is narrowed down to — 1 < to < which in turn implies an interesting 
scenario that universe starts off with a stable state and it retains stability until u> hits 0. Thereafter it will un- 
dergo a stage of thermal instability after which again it will be in a thermally stable state when uj drops below — 1. 

In short, depending on the signature of U\ + 3u;(l + uj) the instability may occur at to = —1 or at to = 
for which the corresponding z is —0.019 & 52.517 respectively. However, there are many objects which are at a 

3 this is the point we do theoretically mark as the phantom barrier, beyond which universe may turn towards different 
singularities like Big Rip(Caldwell et. al 2003; Nojiri et. al. 2005) etc. 
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redshift of z = 52.517. But we are yet to speculate any kind of impacts of thermal phase transition of universe 
upon them. So based on current knowledge of explaining observational data we can rule out such a transition 
at z = 52.517. Besides, the model is valid for z < 4 (Efstathiou, G. 1999). So we can not comment on the 
second possibility staying within the model. 

Note there is observation constraint on ujq as follows: ujq > — 1. So we cant have uj < — 1 unless we make z 
slightly negative i.e we are heading towards instability. Also note u\ is assumed to be positive to reconcile with 
the fact that uj should decrease with time. 

The expression for speed of sound (|14[) within this model also shows instability since with ui going below 
— 1, the square of the speed becomes less than 0. This is generally termed as adiabatic instability. (Bean, R 
et. al. 2008a, 2008b). Bean, R. et al. (2008a, 2008b) concludes the models with nontrivial effective coupling 
between dark matter and dark energy can lead to exponential growth of small adiabatic instability which is 
also characterized by negative sound speed squared. As a result even if universe starts with a uniform fluid 
the instability will bring upon an exponential growth of small density perturbation. In analogy we also can 
speculate that it might be possible that at our indicated points small density perturbation will occur opposing 
fundamental phenomenon like the propagation of sound to take place. 

ASSS Model 

Within ASSS model, the expression for cy depicts that the onset of instability happens when to goes below —2, 
which in turn implies the following timescale of onset of instability: Taking allowed values of parameter (Alam, 
U. et. al. 2004a, 2004b) i.e taking Aq — 0, A\ = — 1.169,^2 = 1.67 we get that the instability happens to set 
on at z = —0.02. The explanation of this scenario will follow the previous discussion. 

6 On the Validity of Thermodynamical Laws 

Study of thermodynamical laws with the universe as a thermodynamical system has been done in any litera- 
ture(Setare, M.R. 2006, 2007a; Setare, M.R., Shafei, S. 2006; Setare, M.R., Vagenas, E. C. 2008; Mazumder, 
N., Chakraborty, S. 2009, 2010; Bhattacharya, S. Debnath, U. 2011). We will just recall the result relevant for 
our models. 

6.1 Validity on Apparent Horizon 

The Friedman metric of a isotropic spatially homogeneous universe is given by: 



From this metric we can easily calculate the radius of the apparent horizon (r' AH ) which comes out to be : 




(63) 



where r' = ar. 



1 



(64) 





Now, in spatially flat universe, k = and thereby we have r' AH = -jj. 

Following Hawking's idea, the temperature associated with the apparent horizon is: 




(65) 



where K is surface gravity of the apparent horizon and given by; K = -r— . 

Hence, the temperature associated with apparent horizon is =—3, — . 

AH 

With k=0, using equation (|7) we obtain 




(66) 
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To determine the validity of first law thermodynamics we first calculate the energy crossing over this apparent 
horizon in an infinitesimal dt time which comes out to be (Bousso, R. 2006) 

4ir(P + p) 3 (l+w) , lR7 , 

dE AH = —^dt = — 2lmT) dz (67) 



On the apparent horizon, the entropy is: 
Hence the following holds : 



Sah = j = (68) 
T AH dS A H = dr' AH (69) 

=*> -dE AH = T AH dS AH (71) 
Hence the first law of Thermodynamics holds on apparent horizon independent of how we are gonna model 

Now, Using Gibbs equation we get (Izquirdo, G. and Pavon D. 2006) 

TdSi = dEi + PdV (72) 
where Si, Ei are entropy and enrgy density respectively inside the apparent horizon and Ei is given by 

E l = ^r AH p = \r AiI (73) 
and the volume bounded by the apparent horizon is 

V = \*r% (74) 

Therefore using equation (|72")l . (j?3"j) and (|75|). we get 

dSi = Trr AH (l + 3oj)dr' AH (75) 

Rewriting the last equation and using the expression for dS A H,~we get 

d(Si + S AH ) 9tt {l+ujf 

dz 2H 2 1 + z [ ' 

^m^a)_^ +uf(x+ , )>a (77) 

Hence holds the Generalized Second Law of Thermodynamics if we consider apparent horizon and the volume 
bounded by it. So the result holds not only for two types of fluid we have discussed in the paper, but also for 
every fluid obeying p = uj(z)p which is in quite agreement with (Mazumdcr, N. et. al. 2010). 

6.2 Invalidity of Thermodynamic laws on Event Horizon 

The event horizon teh is defined to be 

dt 1 f- 1 -1 



r EH = / - = — — / — dz (78) 
a 1 + z J. H 



Hence the event horizon has a hawking's temperature of = — - — which yields 



T E HdS E H = dr EH = \ dz (79) 
(1 + z)H 
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It follows that: 

dE EH + T EH dS EH = 7z „ ( | (1 + w) - 1 ] ^ unless U) = -\ (80) 
(1 + z) i± \l J 6 

Note lu — — i is the marginal point of Strong Energy Condition which is 3p + p > 0. 
Employing similar technique we would get, (Xing, L. et. al. 2011) 

d -^i±^l = 2^ ff f + „ [ 3r |^ 2 (!+») + ?TT ^ (81) 

Changing the variable to a we get; 

d(gj + Sgg) = 2. (r| ff F) (^) + n [3r|^ 2 (1 + w ) + 2r BB ] ± (82) 

Evidently d ^ Sl ^ EH ^ is not necessarily positive for any allowed value of z ( as we have 4^ in the expression) 
which, in turn, implies that the second law of thermodynamics breaks down inside the event horizon in quite 
agreement with Xing, L. et. al. 2011). 

7 Brief Summary 

So far we have studied the thermodynamics of universe from two different models. It has been revealed that 
the first model where lo(z) = ujq + tjjxln(l + z) does not go through any critical point. In this model Universe 
started with a high (fig. 1, 2, 3) energy density, pressure and temperature which falls as z becomes 2 from O.We 
imposed restriction on allowed values of loq and u>\ by demanding v s < 1. The thermal capacity at constant 
volume changes sign as uj becomes smaller than — 1, but this transition is smooth without having critical point. 
Unlike MCG or GCG here pressure depends both on temperature and volume in thermal EoS. From equation 
(|23|) we have shown it is possible to have cy > irrespective of model to ascertain a stable universe. 

Thereafter we have generalized the results for a general EoS parameter. We have shown that energy density 
takes the same form as a function of temperature and volume irrespective of model with a general limiting 
behavior of thermodynamic parameter (|59[) to satisfy third law of thermodynamics. We have also interpreted 
Uq as a measure of number of micro states and argued for it natural plausibility (|6ip and (|62[) . Moreover, it is 
intriguing to note that the thermal EoS (|57|) looks like ideal gas EoS in some way. Also we have determined the 
asymptotic behavior of F(V) such that third law of thermodynamics holds. It has been shown that the first and 
second law of thermodynamics hold on apparent horizon not only for our model but also for any model having 
p = uj(z)p. They fail to be valid on event horizon unless u> takes a specific value which is in perfect agreement 
with (Mazumder, N. et. al.; Xing, L. et. al. 2011). 

One important aspect which have been followed in this paper is a try to speculate the instabilities. Con- 
sidering, the thermodynamical point at which the heat capacity changes its sign, as a transition point we have 
followed that the point is nothing but the phantom barrier in the case of Log parametrization. Even there we 
can find the squared sound speed turns out to be negative. Even for ASSS such point arises for negative z, 
theoretically indicating to some future point. We have also explored the possibility of the fact that at our indi- 
cated points small density perturbation will occur such that the fundamental phenomenon like the propagation 
of sound does not take place at all. 
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